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An iterative random procedure is considered allowing an entanglement purification of a class of 
multi-mode quantum states. In certain cases, a complete purification may be achieved using only 
a single signal state preparation. A physical implementation based on beam splitter arrays and 
non-linear elements is suggested. The influence of loss is analyzed in the example of a purification 
of entangled N-mode coherent states. 
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I. INTRODUCTION 

The impossibility to isolate a physical system com- 
pletely from its environment represents a major diffi- 
culty in the experimental realization of devices process- 
ing quantum information. Of crucial interest are there- 
fore schemes which tolerate a certain amount of envi- 
ronmental coupling yet protect the actual information, 
e.g., by accepting some redundance in the quantum state 
preparations. For example, a pure entangled state may 
be recovered from a number of mixed entangled states 
of spatially extended signals shared by separated par- 
ties A potential application of such entanglement 
purification schemes is the realization of a secure data 
transmission despite using a lossy channel and appara- 
tus 2]. Apart from this, the treatment of entanglement 
purification helps to deepen the understanding of entan- 
glement itself 0, yj . Related questions which have been 
addressed are whether it is sufficient to process each sig- 
nal state copy individually 0] or how to deal with copies 
which cannot be factorized into separate states En- 
tanglement swapping 0,0 has been discussed to propose 
schemes operating more efficient, despite loss. 

While for certain classes of states such a superpositions 
of coherent states, methods solely based on linear optical 
elements like beam splitters and photodetections could 
be found an implementation covering other classes 
of entangled states remains a challenge. Schemes based 
on cross-Kerr couplers have been described for Gaussian 
continuous variable entangled states [9|, llOl lllj . Alterna- 
tively to the generalization to continuous variable states 
shared by two stations, multi-user schemes have been 
considered which involve more than two parties [l^ . 

The aim of this work is to investigate a range of trun- 
cated mixtures of entangled multi-mode states, whose 
periodic nature allows a relatively simple separation of 
their pure components despite the fact that the states 
themselves are infinite. Here, focus is on the theoretical 
possibility of a purification. From a practical point of 
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view, an implementation of the suggested scheme may 
be realistic only for special signal states. The effect of 
loss is discussed in a selected example of application. 



II. INPUT STATES 

Within this work, we assume that the mixture describ- 
ing the input signal state g can be truncated at a suffi- 
ciently large number M, 

M 

0=J2Pn\^n)(^n\, | *n) = <W (1) 

n=0 

By entanglement purification we understand the extrac- 
tion of some pure component |\f f o)( 1 I r o| of g by effectively 
implementing a transformation g' =p~ 1 Y gY^ , p= (Y'Y) 

with Y ~ |4 r )(*ol- It may therefore be understood as a 
multi-mode state purification. We assume that in each 
signal mode j, whose entity is denoted by 1=1, . . . , N, a 
repeater is located whose operation involves some auxil- 
iary modes bj . If the repeaters can only use local prepara- 
tions \F)bj, unitary transformations Uj^. , and detections 
b j (Gj | , they can solely implement local operators Yj (Gj ) 

= Uj(Gj) bj (Gj\Uj,bj \F)bj , which may depend on the re- 
sults Gj of the local detections in the auxiliary modes. 
In contrast, the \^ n ) may be arbitrary entangled states 
of the spatially separated signal modes. Therefore, the 
repeaters must share an entangled state in general. If 
this state has not been provided previously by a separate 
resource similar to teleportation setups, then the only al- 
ternative is to equip each of the repeaters with a second 
signal input mode — j, whose entity is denoted by II = 
— 1, . . . , —N, and which is fed with a second copy of the 
signal state g. If the repeaters periodically receive prepa- 
rations g, the purification can be implemented in a loop 
(using, e.g., a cavity) as a repeated transformation of the 
signal, each time applying a further copy of the signal 
state. 

In what follows, we assume that there exists a single- 
mode basis { allowing the representation of the sig- 
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nal eigenstates in the form of 

oo 

|* n ) = (Af+l)^£q(n)X;W 



1=0 



x \$h(M+l)+ ni )l ■ ■ ■ \®l N (M+l)+n N )N, (2) 

where l = h, . . . ,In, and to ensure normalization we must 
have J2i=o l c '( n )| 2 = 1- The expression J^n^ denotes 



the sum over all (M + 1) N 1 different TV-tuples n = 
m, ■ ■ ■ , Un of numbers rij that can each take the values 



0, . . . , M but must obey [X)j=i n i\ = n - Here, we have 
used the notation 



[m] = Mod(ra, M + 1) 



(3) 



In order to explain the operation of the purification 
scheme, it will be convenient to distinguish between the 
states 



M 



QII 



n=0 
M 



(4a) 
(4b) 



n=0 



where IZi is the state of the signal currently processed in 
the loop and gu that of the original signal just entering 
the second signal input ports of the repeaters. 



III. IMPLEMENTATION OF THE QUANTUM 
REPEATER 

To describe the transformation of the input states 
Eqs. 0, let us take a closer look at the quantum repeater 
illustrated in Fig. ^a). V is a unitary single-mode oper- 
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FIG. 1: Implementation of an iV-mode entanglement purifi- 
cation. The overall setup shown in (b) consists of quantum 
repeaters located in each of the N signal modes I — 1, . . . , N . 
Their operation requires a second copy of the signal state to 
be fed into the inputs II — —1, ... , —N. A detailed view of a 
single repeater is given in (a). 



\n) = V\§n), i.e., the single-mode basis {|$ n )} is supposed 
to be known. A realization of a desired single-mode oper- 
ator or detection in a truncated space also based on beam 
splitter arrays, zero and single photon detections as well 
as cross-Kerr elements is presented in [13J . (Note that the 

V are not required if attention is limited to the photon 
number basis.) Enclosed between devices implementing 

V and V', an array of 2(M+1) cross-Kerr elements is 
placed, whose operation is described by the operator 



(•5) 



This array couples the signal modes j and —j to M+l 
auxiliary modes 0, . . . , M, themselves coupled by two 
2(M+l)-port beam splitter arrays implementing the op- 
erators U and U* , respectively. The latter are defined by 
the matrix elements 



u k i = {ipkWWi) = 
\<Pk) = &£|o)o-- 



i 



VM + T 



e^MTT, (6a) 
(6b) 



Initially, these auxiliary modes are prepared in the state 
\tfo), i- e -i a single-excited Fock state is fed into input port 
of the U -array, while the remaining input ports are left 
in the vacuum. Assuming a final detection of the auxil- 
iary modes in the state | i.e., the photodetectors at 

the output ports m of the W -array detect 6 m ,kj photons, 
the action of the jth repeater without the detector placed 
in the signal output port — j can be formally described 
by an operator 



{<p kj \tfkU\ip ) 



Vf^lU^KUl^V-jVi, (7a) 
1 M 

-r— y 

M + 1 ^ 

m=0 

x fJ-mrtejZj - - k ?) m . (7b) 



The detector placed in the signal output port —j is as- 
sumed to perform a Pegg-Barnett phase measurement 
such that each trial gives a value rrij that can take the 
values 0, . . . , M and is described by the Positive Operator 
Valued Measure (POVM) 



n -j(™j) = ^2\l,mj){l,', 



1=0 



A I 



x \l(M+l) + n)-j. 
In turn, the value rrij determines the operator 



yn = t/.t iTu+rmjQ^ 

J 3 



(8a) 



(8b) 



(9) 



ator that transforms the basis states |<j>„) to Fock states in Eq. {7aJ. 
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IV. INPUT-OUTPUT RELATIONS 

We now consider the overall setup including all signal 
modes as depicted in Fig.^b). If states \tfki) an d POVM 
values m,j are detected in the auxiliary and second sig- 
nal output modes of the quantum repeaters, the reduced 
state of the first signal output modes becomes 



7tj(m,fe) 



Tr„ 



Y(k)Ki <g> gjjyt(fe)n( m ) 



p(m, k) 



(M + 1) 



1-2JV 



p(m, k 

In the first line of Eq. ljTfj|) , k = ki, 



il I W' ik Q I W k . 
kpf and 



N 



Y(k) 



3=1 



(10) 



(11) 



denotes the product of the individual operators Eq. ((Taj- 
Similarly, m = m\, . . . , mjv and 



AT 



n(m) = JJn- J -(m i ) 

3=1 



(12) 



is the product of the respective single-mode POVM's 
Eq. I jSajl . In the second line of Eq. IjlOl) . an operator 
has been introduced to express a rotation of the eigen- 
state indices according to W"t|^ n ) = |*p, l+1 ]), and k — 

\YljLi kj \ . Qi is the original input state gu given by 
Eq. I|4b() but transferred from II to J, i.e., using a for- 
mal mode transfer operator Iiji 

\n) jk (n\ 



Ylf=i h .■• where 



ijk 



n=0 



k (0\e^\0) 3 



P 



(13) 



we may write Qj = I T uQuIj //• Since only the signal 
modes J appear in the second line of Eq. (|10|) , let us leave 
out this mode index in what follows. From Eq. 1)100. we 
see that the state 7Z'(m, k) only depends on k, so that the 
signal output state in case of an event k can be written 
as 



n'(k) 



i 



M 

\ - (fe) 



p{k) ^ ^ 



p(m, k)7Z'(m, k) 



m— k 



p{k) 



(14) 



A repeated run of the purification can then be described 
as an iteration such that the jth iteration step generates 
a signal state 



p(k) 



(15) 



from its precursor originally descending from a 

signal state copy, TZ^ — g. At each step, k — k^ takes 



some value 0, . . . , M. An observer unaware of any detec- 
tion result simply observes the average 



M 



£Cj) = J2p(k)n^{k) =n 



(3-1) 



(16) 



k=0 



so that iZP> = W-^ — g, which is just the unchanged orig- 
inal signal state. Thinking of the complete purification 
procedure as a linear array of repeating stations rather 
than a loop, the resulting setup therefore constitutes a 
'quantum state guide' which only may increase the ob- 
server's knowledge rather than changing the signal state. 
In the special case of an equipartition, p n = (M +1) _ , no 
additional knowledge can be gained even from the mea- 
surements, so that a purification is not possible, TZ^' (k) 
= 72.(3-1), gy taking into account all detection results 
k^\ . . . , few obtained during the steps 1, . the ex- 
plicit expression of the iterated state becomes 



-m fcW,...,^ 



(1) ^TirfcCO 



P [kW,...,kU)] 

M 

52pW\#n){*n\ (17) 



n=0 



with 



P U) 



p [k 



(i) 



-Pn 



-feCO] 



p 



[fc(D, 



i=i 

M 

-Pn\\p 

1=0 



[n-t] ' 



(18) 



where si is the number of events for which k = I. From 
Eq. (fTHfl we see that 72, w only depends on s = sq, ... ,sm, 
so that we can refer the iterated state to this event, 



1 



P(s) 51 p 



k {1) . . k (3) 



p(s) \s 



1 



M 



A I 



n=0 i=0 



p(s) \s 



g s °+ x W^ Q S1 W^ q S2 W^ ■ ■ ■ W^g SM W ] . (19) 



Here ~J2* denotes the sum over all j-tuples k^\ . . . , k^ 
containing s; times the numbers I (= 0, . . . , M), whose 
number is given by the polynomial coefficient P) . Since 
a pure signal input state results in a pure signal output 
state, the random walk of the state may eventually end 
in a pure state. If this happens, the probability of ob- 
taining \^> n ) is p n — (^nlpl^n) as becomes plausible from 
Eq. JT5J}. 
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BINARY OPERATION 



Let us return to Eq. (|15|) . Since we are only interested 
in obtaining the state |^o)i we may limit attention to 
a binary detection by only discriminating between the 
event A; = and its complement fc^O. Eq. i|15fl then gives 
the transformation 



p{0) 



1 



l-p(O) 



(20a) 
(20b) 



Analogous to Eq. 119|) . the explicit expression of the it- 
erated state 



1 



p{q) 



(i- Q y 



(21) 



can be expressed in terms of the number q of events k = 
0. To give an example, in a situation where our desired 
state | "Jo) is still dominant over the impurities, po>p n y 
n^O, the state l^o) is obtained in situations when j~ 1 q 
is sufficiently large. This case is illustrated in Fig. [21 
monitoring on the basis of the expectation value of the 
von Neumann entropy S = — In TZ^ a successful exam- 
ple of a simulated random purification process calculated 
according to Eqs. H20|) . 
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FIG. 2: Typical evolution of the von Neumann entropy S = 
-Tr[ft (j) ln"fc (i) ] of the signal state ti u) as seen by an ob- 
server aware of the measurement results occuring during the 
first j cycles of a left alone (random) purification process ac- 
cording to Eqs. 12011 . The initial probabilities are assumed to 
be p n ~ e _n , where n = 0, . . . , 100. 



VI. PURIFICATION IN A SINGLE INSTANT 

If the signal state is a mixture of the (M + 1) N states 



|*n) = ^d(n)|$, l(M+1)+ni }i JV (M+l)+njv }n 

(22) 



1=0 



instead of Eq. J5J , each copy can be purified individually 
without need of further inputs. All we have to do is feed- 
ing the second signal input ports of the repeaters with 
vacuum states |0)_j and removing the devices imple- 
menting V-j thus making the -measurements redun- 
dant. Together with £j = En=oP™l*™><*™|, E T is 
then replaced with 



e'Ak) = 



Ti- 



ll 



Y(k)Q! (8 |0) ZJ <0|ft(fe) 



p(k) 



= |*fc)<* fc |. (23) 
The probability of obtaining |$fe) is just pk. 



VII. 



INFLUENCE OF LOSS IN AN EXAMPLE 



Let us return to mixtures of states Eq. (J2J. From a 
practical point of view, the loss occuring in the purifica- 
tion loop must be taken into account. The problem is 
that its effect depends on the respective circulating sig- 
nal state TZ itself. Let us therefore choose the particular 
example of a superposition 



g(r) 



(\a)(a\ + \-a)(-a\) +r(\a)(-a\ + \-a)(a\) 



2(l + re -2|«l 2 ) 



i 

J2p n \^n)(^r, 
n=Q 



(24) 



of two iV-mode coherent states \a) — \oti)i • • • \cxn}n, 
characterized by their complex amplitude a=ai, . . . , ajv 

according to \ctj)j =e aja j~ a i aj \0)j. The eigenvalues and 
eigenstates are given by 

1 + C-iy i e" 2|a|2 
Pn = 1\ „!_,,^ + (25a) 



2(1 + rc- 2 l a l : , 
jo) + (-l)"| -a) 

N 



(25b) 



where we have defined |a| = X)j=i \ a j\ ■ ^ n the meso- 
scopic case, e 
form of Eq. (J2J), 



1*0 



3 = 

2|a, | ^ ^he jij,^ can written in the 
2^ V (n) l$„ 



|q 3 -) + (-!)"' !-<*,-) 
v/2[l + (-l)^e- 2 l^l 2 ]' 



(26a) 
(26b) 



The pure states l^n), to which Eq. (|24|) reduces for r 
= (—1)™, just represent the superpositions |0, ...,0) ± 
|1, . . . , 1) considered in [l2|. Modeling the decoherence 
of the state Eq. (|2*4*jl by the master equation 



L 



do 

dx 



N 

Vj ( 2 «j Qa] -fijg- ghj ) , 

i=i 



(27) 
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where x may be, e.g., the propagation distance, we see 
that the purity parameter r and the individual complex 
amplitudes otj decrease according to 



r(x) = r(0)e 2 ( 



a| 2 -|a | 



2\ Z«L 



r (0)e- 4 l Q °l (28a) 



a (x) = aj(0)e ,,J ^ 



«j(0), 



(28b) 



where r\— X^^li Vj i s determined by the rjj > defining 
the (e.g., scattering) losses and ao = a(0). For distances 
small compared to the classical transparency length, x 
<C L, the damping of the complex amplitudes can be ne- 
glected. For given complex amplitudes, the state Eq. 124(1 
is then determined by its purity parameter r. 

The aim of our purification process is therefore to re- 
obtain a state with r = 1 from a number of states with 
< \r\ < 1. Our model of a lossy purification is that dur- 
ing the feedback between subsequent purification steps, 
the signal state undergoes a decoherence according to 
Eq. l|27jl. which leads to a damping of the purity pa- 
rameter by some factor rfw < 1, which we call feedback 
efficiency. To obtain the resulting evolution of r, we in- 
sert Eq. into Eqs. i.e., H^-^=q[R^-^] and g 
— 8( r ), which gives a state g[R^} given by Eq. I|24|) with 



If 



1 iri?^- 1 ) 



(29a) 
(29b) 



Here, a positive sign corresponds to the event '0' and 
negative sign to the event £ -i0' . With the stationary 
condition flW = R^~^ = R (oo \ Eqs. Q give a lower 
and upper bound 



) - ?]f) 2 + 4r 2 7yF - (1 - t?f) 



2r 



(30) 



between is able to walk. To allow an increase of 

|i£W| above its initial value |r|, the condition > 
|r| must hold, which gives 



1 + 



(31) 



It follows that the required feedback efficiency increases 
with the initial purity. For r/p < 1/2, the scheme is of no 
use at all. On the other hand, the ideal value R(°°) = 1 
can only be approached in practice. In order to keep the 
bound sufficiently close to one, |i?(°°)|>l-e with some 
given £<C 1, the requirement 



riF I 1- 



2e\r\ 
1 + Ir 



(32) 



must hold. For sufficiently small e, a behavior similar to 
the perfect case is observed. Note that the case consid- 
ered here is analogous to the example of a related two- 
mode state discussed in |l4|. For a more detailed analy- 
sis including numerical simulations, we therefore refer to 
that work. 

VIII. CONCLUSION AND OUTLOOK 

We have described a procedure allowing a purification 
of a class of iV-mode quantum states as an iterative ran- 
dom process. While in general, a number of identical 
signal state preparations is applied, only a single prepa- 
ration is required in certain cases of mixed states. The 
physical implementation suggested involves beam split- 
ter arrays, zero and single photon detections as well as 
cross-Kerr elements. The role of imperfections of the 
procedure is modeled in the example of a superposition 
of two TV-mode coherent states. In this case, the de- 
mands on the device increase with the initial purity of 
the input state to be enhanced. A main drawback of the 
scheme is its relying on large non-linearities which are 
difficult to implement at present. An alternative solely 
based on, e.g., beam splitter arrays and photodetections 
would therefore be desirable. From a theoretical point of 
view, the scheme suggests that different types of entan- 
glement may also be distinguished by the effort required 
for a purification. 
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